Exchange Boson Dynamics in Cuprates: Optical Conductivity of HgBa 2 Cu04 + <5 
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The electron-boson spectral density function I 2 x(&) responsible for carrier scattering of the high 
temperature superconductor HgBa2Cu04+a (T c — 90 K) is calculated from new data on the optical 
scattering rate. A maximum entropy technique is used. Published data on HgBa2Ca2Cu308+i (T c = 
130 K) are also inverted and these new results are put in the context of other known cases. All spectra 
(with two notable exceptions) show a peak at an energy (O r ) proportional to the superconducting 
transition temperature Q r ~ 6.3I;bT c . This charge channel relationship follows closely the magnetic 
resonance seen by polarized neutron scattering, Q" eutron m 5.4fesT c . The amplitudes of both peaks 
decrease strongly with increasing temperature. In some cases, the peak at Q r is weak and the 
spectrum can have additional maxima and a background extending up to several hundred meV. 

PACS numbers: 74.25.Gz, 74.62.Dh, 74.72.Jt 



Superconductivity arises when electronic quasiparti- 
cles bind together into Cooper pairs which condense into 
a macroscopic quantum state. The interaction which 
drives the pairing could be charge or spin polarizations, 
but in conventional metals, it is the polarization of the 
lattice of ions which provides the necessary attraction. 
There is also a smaller Coulomb repulsion. In the phonon 
case Eliashberg theory accounts for all the experimen- 
tal data with its central ingredient, the electron- phonon 
spectral density a 2 FjXl) specifying the boson exchange 
pairing interaction Optical spectroscopy is a pow- 
erful tool that can be used to extract bosonic spectral 
density with high degree of accuracy in a variety of sys- 
tems where other techniques such as tunneling and angle 
resolved photoemission (ARPES) arc difficult to apply. 
It has been used in conventional metals @, 0] and in the 
high T c oxides 0, @, @, 0, S, H, E3] but its interpretation 
in the cuprates remains controversial. As emphasized 
recently by Anderson [ill ], the pairing interaction could 
involve high energy virtual transitions across the Mott 
gap with energy set by the Hubbard U which is a large 
energy and the effective interaction would be instanta- 
neous on the time scale of interest. On the other hand, 
recent numerical work 12, [l3[ based on the t-J model, 



has shown that the main contribution to the pairing glue 
is provided by the spin fluctuations with characteristic 
energies of at most a few hundred meV. Optical experi- 
ments provide a direct probe of this energy region. 

The mercury versions of the cuprates, HgBa 2 Cu04 +( 5 
(Hgl201) and HgBa 2 Ca 2 Cu30 8+( 5 (Hgl223) provide a 
unique opportunity to test these ideas. The one-copper 
layer system Hgl201 shares a conventional transition 
temperature around 90 K with widely studied systems 
YBCO and Bi2212, whereas the three-layer compound 
Hgl223 has a T c = 130 K, nearly 40 % higher. These 



dramatically different T c 's lead, as we will show, to very 
different bosonic spectra and place severe constraints on 
models of superconductivity in the cuprates. 

In this letter, we present new data on the optical con- 
stants for HgBa 2 Cu04 +( 5 (Hgl201). The bosonic spectral 
density J 2 x(f2), recovered by maximum entropy inver- 
sion, is found to have a remarkable resemblance to pre- 
vious results for optimally doped Bi2212 [8|. For further 
comparison, using the same methods, we also invert pub- 
lished [IH optical constants in the three-layer Hgl223. 
The high-quality Hgl201 single crystal was grown by a 
flux growth technique. Our sample was a millimeter sized 
platelet with a well oriented a6-plane. It is slightly un- 
derdoped with a T c = 91 K. The real and imaginary parts 
of the optical conductivity a(T,tu) = <?i(T, oj) + icr 2 (r, u>) 
follow from the reflectivity. In analyzing optical data, we 
use a memory function or optical self energy E op (T, u>) 
instead of working with <r(T. k>)[l5l Il6j |. By definition, 
a(T,uj) = (iw 2 /47r)/[tj - 2£°p(T,w)] where uj p is the 
plasma frequency. The optical self energy defined this 
way plays a role analogous in optics (which involves 
a two particle process) to the quasiparticle self energy 
Y, qp {T,Lu) in ARPES [Hill- The optical scattering rate 
is 1/t°p{T,uj) = -2E 2 p (T,w) and the optical effective 
mass m op {T,Lu)/m (with m the bare electron mass) is 
given by uj[m°P(T,uj)/m - 1] = -2E° P (T, uj). 

Our results for the real part of the optical self en- 
ergy — 2E° P (T, uS), which provides an easy comparison 
with ARPES, often presented in terms of 2Hf(T, cu), are 
shown in Fig. [T]as a function of ui for eight temperatures. 
The lowest three temperatures are in the superconduct- 
ing state. The most prominent feature of the curves 
is the peak around 100 meV seen at T = 29 K, which 
progressively disappears as the temperature is increased. 
To better see how the bosonic spectrum generates the 
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FIG. 1: (color online). Experimental results for the real part 
of the optical self energy — 2£° P (T, u>) as a function of pho- 
ton energy u for eight temperatures. Inset, theoretical re- 
sults based on numerical solutions of the generalized Eliash- 
berg equations. (Solid blue superconducting and dashed black 
normal state.) The electron-boson exchange spectral density 
used is shown as the dashed-dotted blue curve. The super- 
conducting gap value is A = 22.4 meV. 



self energy we show in the inset of the figure a model 
calculation based on numerical solutions of the d-wave 
Eliashberg equations. The input electron-boson spec- 
tral density I 2 x(Q), shown as the dashed-dotted curve, 
consists of a large narrow peak (right hand scale) cen- 
tered at 56 meV followed by a dip and a long back- 
ground extending beyond 400 meV. The dashed curve, 
with the same bosonic spectrum, is the normal state re- 
sult for — 2£° P (T, u>) which is to be compared with the 
solid curve in the superconducting state. The supercon- 
ducting gap obtained was 22.4 meV giving a gap to T c 
ratio, 2A/fcsT c = 5.8. Note that, as theory predicts [l9j |. 
the dashed curve in the normal state shows no visible 
structure at u> = f2 r = 56 meV. Instead there should be 
zero slope at uj = y2fl r for an Einstein spectrum. It is 
clear that boson structure is hard to see in normal state. 
However in the superconducting states the quasi-particle 
electronic density of states acquires energy dependence 
and this helps reveal the underlying boson structure as 
seen in the solid blue curve. The peak at ~100 meV is 
neither at Q r nor \/2D, r but is shifted upwards by the 
opening of the gap A [l^ but as we see, the position of 
the peak in I 2 \(^) cannot be read off the curve without 
the knowledge of the value of the superconducting gap. 

In Fig. [DJb) we show results for the spectral density 
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FIG. 2: (color online). Top frame, the optical scattering 
rate 1/t° p (T,u>) for Hgl201 v.s. uj for 8 temperatures (light 
curves). The wider curves are our maximum entropy recon- 
structions. Bottom frame, the electron-boson spectral func- 
tion I 2 x(ty v.s. Q. The inset gives the peak position (blue 
triangles) left scales as a function of temperature and the red 
squares give the corresponding peak amplitude. 



7 2 x(J7) of maximum entropy inversions augmented with 
a least squares improvement based on the full d-wave 
Eliashberg equations [6J]. Further applications are found 
in Refs. [1, The input to the inversion is the op- 

tical scattering rate 1/t{uj) — -^-TZe(l/a(uj)). These are 
shown in Fig. a) where the results of our inversions 
(heavy lines) are compared with the original data (light 
lines). In all cases, the fit is very good. Recently, van 
Heumen et al. 2^, 21 1 have also presented optical data for 
Hgl201 above T c which they analyze in terms of a spec- 
tral density represented by a set of histograms. While 
they obtain fits which are of equal quality to ours and 
have a background extending to > 400 meV as we have, 
they find that the height of the peak at 56 meV does not 
change contrary to our findings for T = 246 and 295 K. 
This has been taken as evidence for coupling of the charge 
carriers to phonons [22|. We can also get fits where the 
peak height does not change but only if we use a biased 
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FIG. 3: (color online). 
(T c =130 K). 



Same as Fig. [2] but for Hgl223 



maximum entropy inversion with the default model set to 
the previous lower temperature solution instead of being 
set to the constant of the unbiased inversion. 

In the inset to the lower frame of Fig. [2) we show the 
frequency (left scale, triangles) of the prominent peak in 
I 2 x(ty as a function of temperature. As also noted by 
van Heumen et al. [20], Q r is fairly constant at ~ 56 
meV but in our analysis, this frequency clearly increases 
for T above 200 K. More importantly, the amplitude of 
the peak shows strong temperature dependence in the su- 
perconducting state and also above 200 K. Further, the 
width of thepeak increases with increasing T. As noted 
for Bi2212 [8(, the shift in spectral weight into the peak 
at 60 meV can be interpreted to proceed through a trans- 
fer of spectral weight from high to low frequencies as the 
temperatures is lowered. We note here that in contrast 
a bosonic function from the electron phonon interaction 
would not have these properties: its amplitude, width 
and center frequency would all be temperature indepen- 
dent and would not vary from one cuprate to another (as 
shown in Fig. 4 below). 

In Fig. [3] we display similar results for Hgl223 with 
T c = 130 K. In Fig.^a) we reproduce the optical scatter- 
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FIG. 4: (color online). The optical resonance frequency Q r 
as a function of T c . Bi2212 @,i], T12201 H, YBCO [llflQ, 
LSCO 9], and PCCO [T(|. 



ing rate at four temperatures (with T = 15 K and 125 K 
in the superconducting state) from the work of McGuire 
et al. [14]. Our results for I 2 x(Q) are presented in Fig. 
[3f b) and the quality of the data reconstruction is demon- 
strated by the heavy lines in Fig. El^a) to be compared 
with the corresponding light line (data). We find signifi- 
cant residual (static impurity 1 /rj mp ~ 95 meV) scatter- 
ing rate in contrast to Hgl201 which is in the clean limit 
(l/ T im P = 0). The superconducting state data (blue, 
dotted curve) for l/r op (w) at T = 15 K shows a peak 
around 140 meV which is the indication of a gap in the 
charge carrier density of states (DOS). Currently there is 
no known kernel which allows maximum entropy inver- 
sion of such data and so we do not show results in this 
instance. On the other hand, in the T — 125 K data there 
is no signature of such a DOS-gap. Results are shown in 
Fig. [3] b). The prominent peak at fl r = 72 meV seen 
in the curve at T = 125 K is missing at higher T. In 
contrast to the HgT201 case no reasonable alternate fits 
can be found for T = 225 K and 295 K which show a 
significant peak amplitude at 72 meV. 

In Fig. [U we place our results for f2 r , the frequency of 
the peak in I 2 \(Q), in the context of other such results 
by plotting f2 r as a function of the superconducting T c 
for a number of cuprates. In all cases, I 2 x(f2) has been 
extracted from the optical data but not always using a 
maximum entropy technique. Some are fits to assumed 
forms including a broad background introduced in refer- 
ence [23} to model antiferromagnetic fluctuations, aug- 
mented with a resonance peak. Both methods give very 
much the same results as documented in reference [5j. 
The heavy long dashed line is a least square fit to all the 
optical data and gives fl r « 6.3fcsT c . This is close to, 
but not quite, the position of the spin-one neutron res- 
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onance obtained by He et al. [24], [2J| where n™ utron w 
5.4fcsT c represented in Fig. [4] as the dotted line. 

Several comments should be made about such a com- 
parison between charge excitations and the magnetic sus- 
ceptibility. First, the the neutron resonance plotted in 
Fig. prefers to the sharp peak that appears at q = (ir, n) 
whereas the bosonic function that governs the optical re- 
sponse is the q averaged local susceptibility. This arises 
because the Fermi-surface to Fermi-surface electron scat- 
tering involves momentum transfers to boson excitations 
that span all momenta in the Brillouin zone, some involv- 
ing Umklapp processes. Where magnetic neutron scat- 
tering data for the g-averaged susceptibility are available 
such as the Ortho II YBCO [2(| or optimally doped 
LSCO 27] the agreement between the neutron data and 
the optical data is excellent: not only are the peaks in 
the response at the same frequencies but also the tem- 
perature dependence of the amplitude of the peaks are 
in agreement [1, 0] ■ Secondly, it should be noted that in 
some cuprates, the resonance described involves only a 
very small fraction of the total weight seen in the local 
spin susceptibility as in LaSrCuO Q and PrCeCuO [lj| • 
In YBCO6.5 it is estimated to be 3% (26| and in some 
cases there is no resonance 0, l28| . but it is always the 
local, Brillouin zone averaged spin susceptibility which 
controls superconductivity. Finally, the two points at fl r 
= (not used in the fit to the data in Fig. [3J are for 
YBC0 6 .35 (HI and overdoped T12201 @. In both cases, 
no optical resonance could be identified. The resonance 
may enhance but is not essential for superconductivity 
in the cuprates and the scaling of the position of the 
peak with T c shown in Fig. 0]must be the result and not 
the cause of the rearrangement of the electronic DOS in 
the superconducting state as suggested by several theo- 
rists [3, [3(| ■ We also note here that recent dynamical 
mean field calculations of the one-band Hubbard model 
yield bosonic spectral functions very similar to what is 
shown in Figs. 2 and 3. 12, 13| 



In summary we find that in Hgl201 and Hgl223 opti- 
cal resonances are found in maximum entropy inversions 
of the optical scattering, at 56 and 72 meV, respectively. 
However, when the temperature is increased towards 300 
K, the spectral weight under this resonance moves to 
higher energy and broadens significantly, in contrast to 
the findings of van Heumen et al. 21|. The optical res- 



onance scales with T c over a broad set of materials with 
£! r w 6.3fcsT c which is remarkably close to the energy 
of the spin one resonance seen in polarized neutron scat- 
tering, namely $} neutron = 5.4fc B T c leaving no doubt that 
the charge carriers are coupled to spin fluctuations, while 
there is no evidence for an important phonon contribu- 
tion. 
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Note added in proof. We have learned of a neutron scat- 
tering study 4 by Yu et al. [31], where a magnetic reso- 
nance in optimally doped Hgl201 is reported at 56 meV, 
exactly the same energy as the peak we found here. Re- 
cent Raman data 5 nd a superconducting gap 2?? very 
close to the values found in our calculations [32]. 

Note added in proof. -We have learned of a neutron 
scattering study by Yu et al [31], where a magnetic reso- 
nance in optimally doped Hgl201 is reported at 56 meV, 
exactly the same energy as the peak we found here. Re- 
cent Raman data finds a superconducting gap 2 A very 
close to the values found in our calculations. [32l 
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